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Abstract 

We examine the sensitivity and surface bias of convolved DGLV radiative and elastic 
loss in a realistic geometry including Bjorken expansion. We find that this more 
faithful treatment of the medium density is not a priori reproducible via fixed length 
approximations, and neither the fragility nor the surface emission of BDMPS-based 
models is seen. 
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1 Introduction 

The ultimate goal of hard probe physics is jet tomography, the use of the high- 
p± hadronic attenuation pattern to learn about the medium density. In order 
to believe that this is an achievable goal one needs a theoretical description of 
the processes involved that describes the data and that has its parameter(s) 
relatable to the matter density under calibrated control. 

For these tools to be truly tomographic, however, there is the additional re- 
quirement that the resulting medium measurement be reasonably precise. The 
overall precision of density determination is set by a combination of experimen- 
tal precision and theoretical sensitivity: if the experimental data are infinitely 
precise, any model with the previous paragraph's characteristics could asses 
the medium density with similar precision; but experimental error bands exist 
and must be taken seriously, limiting the precision of medium determination 
by the innate sensitivity of the theoretical model to changes in the afore- 
mentioned parameter (s). As time progresses experimental error bars tend to 
decrease, allowing for a more precise measurement. 
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[1] claimed that BDMPS-based radiative energy loss models are a fragile probe 
of the medium; i. e. large changes in q, the input parameter for the theory, are 
not reflected in Raa predictions incompatible with data. Surface emission, the 
process in which only the partons produced very close to the medium edge 
escape for observation, was posited as the reason for their observed insensi- 
tivity to increased medium opacity. This is an appealing argument, and the 
latter two of [1] provided evidence for this idea. The very pessimistic conclu- 
sion, then, was that jet tomography is not possible at RHIC and will not be 
possible at the LHC. 

We counter the dire claims in [1] by first noting that several approximations 
were made in their description of the medium produced in heavy ion collisions. 
In the first of [1], the authors used a hard cylinder initial nuclear density and 
did not include Bjorken expansion. The latter works of [1] used the realistic 
Woods-Saxon nuclear density distribution but again neglected Bjorken ex- 
pansion. These certainly aided in alleviating some of the high numerical cost 
involved in the calculations. However we feel these geometrical simplifications 
strongly biased their results to the surface of the medium resulting in a signif- 
icant loss in model sensitivity. Moreover, experimental error bars have shrunk 
considerably in the time since [1] were published. These issues raise serious 
doubts about the validity of their conclusions. In this paper we will explore 
the sensitivity and surface bias of a theoretical model based on radiative and 
elastic energy loss in a realistic diffuse medium that includes Bjorken expan- 
sion. 



2 Results 

Unfortunately, the use of realistic geometry is nontrivial. The analytical radia- 
tive energy loss calculations necessarily were made in simpler geometries, and 
the elastic energy loss involves a line integral through the local temperature 
field of the medium; a mapping from the diffuse, p part -like medium into the 
approximate, idealized geometrical settings in which the theoretical calcula- 
tions were made is required. For the WHDG model using convolved elastic and 
inelastic energy loss, we define the parton's path "length" as L inpnt (x,n; b) = 
I dlp paxt (x + n/;6)/(p part )(6), where (p paTt )(b) = J ' dxdyp^ ait (x,y;b) /N pait (b). 
The probability of a parton having a particular path length is then given 
by P(L; b) = J dxdyd(f)TAA(x] b)S(L — L input (x, n; b)/2nN pa _ rt (b). We consider a 
diffuse Woods-Saxon participant density created by the Glauber profiles and 
include ID Bjorken expansion. Fig. 1 (a) shows the wide distribution P(L) 
of input path lengths resulting from a most central collision. We also find 
the flavor-dependent single lengths Lq, Q = g,u,c,b that best reproduce the 
energy loss resulting from the full P(L) distribution; these are represented 
by the vertical lines in Fig. 1 (a). What is most striking in Fig. 1 (a) is the 
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hierarchy of Q-dependent length scales. No single, representative path length 
can account for the full distribution of lengths, and their hierarchy of lengths 
is not a priori obvious. 
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Fig. 1. All plots depict midrapidity for most-central collisions, (a) Distribution of 
path lengths traversed by hard scatterers; a purely geometric quantity, it is the same 
for all jet varieties. Note that no single pathlength can reproduce the results of the 
full distribution for all parton flavors, (b) Model sensitivity to varying dN g /dy , 
which is much greater than observed in [1]. Susceptibility to changes in a s suggests 
the need for a relaxing of the fixed coupling constant approximation, (c) Transverse 
coordinate (x, 0) distribution of surviving p± = 15 GeV, Q = g,u,c,b jets moving in 
direction <j> = as indicated by the arrows. Units are arbitrary for illustration. The 
binary collision distribution of initial jet production points, pj e t(x, 0), is shown; the 
ratio pq I pj e t gives the local quenching factor including elastic and inelastic energy 
loss though the bulk QGP matter distributed as pqqp(x,0). 

A detailed explanation of our DGLV radiative convolved with elastic energy 
formulation can be found in [2]. The single input parameter is dN g /dy, which 
sets the normalization scale between the actual medium density and p pa rt- We 
mention briefly that, as discussed in [2], we have a fixed a s coupling. One 
can see from Fig. 1 that our model agrees with the data for a medium den- 
sity compatible with the observed entropy constraint, dN g /dy ~ 1000. Using 
fluctuating path lengths in a realistic and expanding background medium ge- 
ometry, the predicted Raa is no longer consistent with data when the medium 
density is artificially enhanced by a factor of two. Our sensitivity to medium 
density changes is enhanced over pure radiative loss models by the inclusion of 
elastic energy loss, due to its smaller width for fluctuations relative to radiative 
[2]. We note, however, that models using GLV-type radiative loss or higher- 
twist loss only are similarly not fragile [3]. We therefore conclude that the 
pion Raa(p±) is in fact quite a good experimental signal for jet tomographic 
studies at RHIC. 



We also do not observe surface emission in our model of RHIC physics. De- 
fine a parton's survival probability, the chance that a parton produced at a 
specific location and making an angle with respect to the reaction plane is 
observed escaping the medium, as Pq(x, 0) = pj ot (a ? ) / de(l — e) n Pg(e; L(x, 0)), 
see Fig. 1. Again one sees that no single, representative path length can ac- 
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count for the distribution of all flavors. While all are biased to the surface, 
a necessary consequence of nonzero energy loss, none of the distributions can 
be categorized as surface emission. The characteristic widths of these distri- 
butions range from Ax ~ 3 — 6 fm. Similar results of partonic survival from 
deep within the medium are seen in more careful BDMPS-based studies that 
use realistic geometry and Bjorken expansion and the higher-twist energy loss 
model, [4]. We therefore conclude that the observed pions do explicitly probe 
significant portions of the medium. 



3 Conclusions 



In order to make strong statements on fragility, energy loss models must use 
realistic nuclear densities with Bjorken expansion. Otherwise the naive expla- 
nation of model insensitivity to increased medium density, surface emission, 
is a reflection of the surface bias of the oversimplified base geometry, not of 
the actual surface tendencies of observed high-p^ partons. Unfortunately, the 
use of realistic geometries is nontrivial; RHIC is not a brick. Our model for 
high-p^ attenuation with convolved elastic and inelastic energy loss in a real- 
istic geometry with Bjorken expansion is quite sensitive to changes in medium 
density and does not display the properties of surface emission. 

Further work is needed: as evidenced by Fig. 1, R A a depends strongly on the 
coupling; the assumption of fixed a s must be relaxed. The first reference of [4] 
showed that several trivial unphysical models approximated the nearly flat in 
Raa data; we must show that our model of radiative, elastic, and path length 
fluctuations can reproduce the plethora of high-p^ observables. If we convince 
ourselves that our theoretical description is correct, then we may return to the 
precise experimental Raa measurement to tomographically probe the medium. 
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